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ABSTkACT 


The  wind  field  is  investigated  which  is  encountered  by  a  missile 
traveling  along  a  hypothetical  trajectory  downwind  of  a  two-dimensional 
ridge.  Reasons  are  given  for  studying  this  situation  in  a  wind  tunnel. 

The  problem  is  reduced  to  the  determination  of  turbulence  spectra  and  of 
joint  probabilities  for  the  joint  occurrence  of  two  velocities  simultane¬ 
ously  along  the  trajectory  which  corresponds  to  mean  flow  conditions. 

In  the  theoretical  part  an  attempt  is  made  to  obtain  approximations 
to  the  joint  probability  density  distributions  which  yield  to  experimental 
evaluation.  The  experimental  part  is  concerned  with  measurements  of  pro¬ 
files  of  mean  velocities  and  turbulent  intensities  and  with  the  -'termin¬ 
ation  of  turbulence  data  for  evaluating  spectra  and  joint  probability 
distributions . 


A  preliminary  version  of  this  report  has  been  presented  at  the  Unguided 
Bai’istic  Missile  Meteorology  Conference,  Las  Cruces,  Oct.  ol  -  Nov.  2, 
.  :<b7  . 
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Chapter  I 
INTRODUCTION 


1,1  The  Problem 


One  of  he  major  problems  in  predicting  the  target  hitting  capabilities 
of  unguided  rocket  propelled  missiles  flying  in  the  atmospheric  boundary 
layer  is  the  interaction  between  the  missile  and  the  turbulent  wind  field 
along  its  flight  path.  In  the  analysis  of  missile  weapon  systems,  es¬ 
pecially  those  used  in  short  range  (0-1  km)  applications,  predicting 
target  hit  probability  caused  by  gust  winds,  involves  prior  knowledge 
„ f  the  wind  field  along  the  missile's  trajectory.  We  can  formulate 
tnis  problem  as  follows:  if  the  trajectory  of  a  missile  is  given  by  a 
deterministic  curve  determined  by  mean-wind  conditions,  we  must  find  the 
probability  distribution  of  the  perturbations  of  the  trajectory  end  point 
if  the  missile  encounters  random  velocity  fluctuations  during  its  travel 
along  the  trajectory.  The  fluctuations  influence  the  flight  path  in  two 
ways.  Vibrations,  caused  by  the  gust  spectrum  might  occur,  and  the  missile 
might  be  deflected  from  its  course  by  large  velocity  fluctuations.  For 
obtaining  instantaneous  wind  measurements  to  calculate  trajectories  in  a 
turbulent  wind  field,  the  present  experimental  study  was  undertaken. 

we  chose  the  wind  tield  which  exists  in  the  wake  downwind  of  a  two- 
dimensional  obstruction  with  air  flow  separation  at  the  downwind  slope, 
as  shown  in  Fig  1.  The  sinusoidal  obstruction  used  in  this  study  repre¬ 
sents  the  model  of  a  ridge.  The  wind  field  which  exists  in  the  wake  of 
a  ridge  is  of  interest  in  military  co»ibat  applications  since  ridges  have 
been  used  as  part  of  a  defensive  line  against  an  attacking  force.  If 
missile  launchers  are  emplaced  ah-.,g  a  ridge,  the  target  impact  dis¬ 
persion  of  missiles  caused  by  the  turbulent  winds  on  the  lee  side  of  the 
ridge  will  play  a  considerable  role  in  battlefield  strategy. 

A  full  account  of  this  wind  field  is  difficult  to  obtain  in  the 
field.  The  number  of  data  points  at  which  wind  speed  information  is 
required  is  large,  and  the  variability  of  wind  speeds  in  natural  envi¬ 
ronments  would  require  elaborate  and  costly  experimental  equipment. 
Therefore,  it  was  suggested  to  study  the  wind  fields  that  might  be 
encountered  downwind  of  a  sinusoidally  shaped  hill  in  the  controlled 
environment  of  a  laboratory  where  many  needed  data  can  be  taken  one 
after  another  instead  of  simultaneously,  and  where  the  reliability  of 
measuring  instruments  and  data  analysis  equipment  has  reached  a  high 
level . 


In  this  report,  wc  shall  concentrate  only  on  the  problem  of  ob¬ 
taining  an  approximation  to  the  joint  probability  distribution  for 
a  sequence  of  instantaneous  velocity  vectors  along  some  hypothetical 
trajectories.  The  analytical  considerations  are  based  on  assuming  cer¬ 
tain  models  for  joint  probability  distributions.  The  validity  of  these 
distributions  for  the  disturbed  flow  field  downstream  of  ..  ridge  is 
demonstrated  by  means  of  experimental  data  obtained  in  the  wind  tunnel. 
The  observations  were  made  for  a  steady  mean  velocity  field  obtained  bv 


blowing  air  parallel  to  a  flat  plate  perpendicularly  onto  a  model  ridge 
of  sinusoidal  shape. 

1  .2  Considerations  on  Modeling 


The  crucia:  problem  in  applying  laboratory  results  for  practical 
applications  in  a  natural  environment  is  the  question  of  scaling  labora¬ 
tory  conditions  up  to  field  dimensions.  For  flows  of  undisturbed  bounda¬ 
ry  layers,  such  as  the  wind  along  a  boundary  of  constant  roughness  over 
a  long  fetch,  the  modeling  has  been  achieved  beyond  reasonable  doubt  by 
scaling  according  to  the  ratio  of  the  roughness  heights,  and  by  keeping 
the  shear  velocities  constant.  With  these  conditions  met,  both  the  mean 
velocity  conditions  and  the  turbulence  structure  are  approximately  scaled 
For  a  boundary  layer  flow  which  is  disturbed  by  a  sharp  edged  obstacle, 
Plate  and  Lin  (1965)  have  presented  an  argument,  based  on  the  boundary 
layer  integral  momentum  equation,  that  the  same  parameters  together  with 
the  drag  coefficient  of  the  obstacle  (as  referred  i  '•  some  convenient 
velocity,  such  as  the  geostrophic  wind  velocity),  suffice  to  model  the 
mean  velocity  field.  As  far  as  the  turbulence  structure  is  concerned, 
no  equivalent  conclusions  are  as  yet  forthcoming,  but  some  work  by 
Plate  and  Lin  (1966)  has  pointed  at  the  possibility  that  the  modeling 
of  the  dissipation  number  is  an  additional  requirement.  Moreover,  no 
conclusions  have  yet  been  reached  on  how  the  turbulence  structure  would 
be  affected  if  this  number  is  not  modeled  accurately.  Work  is  in  progres 
on  this  point  at  Colorado  State  University.  It  is  reasonable  to  suspect 
that  modeling  requirements  will  result  in  a  scale  factor  for  the  dissi¬ 
pation  rates  which  does  not  differ  very  much  from  that  for  the  mean 
velocity. 

With  this  assumption  made,  translation  of  laboratory  data  to  field 
data  is  a  simple  problem,  provided  that  the  drag  coeificient  of  the 
obstruction  can  be  estimated.  The  procedure  would  be  to  determine  the 
roughness  length  and  the  geometrical  pattern  of  the  natural  situation, 
and  then  to  prepare  a  scale  moo  1  of  it  in  the  laboratory,  setting  the 
roughness  length  in  the  laboratory  at  a  convenient  level  by  artificial 
roughening  of  ihe  wind  turn  el  boundary  As  long  as  the  dimensions  of 
the  obstruction  are  such  that  it  lies  well  within  the  lowest  1000  to 
2000  ft  of  the  atmosphere,  and  as  long  as  the  wind  velocity  is  such 
that  the  gross  Richardson  number  of  the  prototype  is  not  essentially 
different  from  zero,  and  as  long  as  the  model  is  sharp  edged,  so  that 
the  separation  line  is  fixed,  the  condition  in  the  laboratory  should  be 
similar  to  that  in  the  field 


°model  °field 

In  this  equation,  h  is  the  height  of  the  obstruction  and  z  is  the 
roughness  height. 

For  an  obstacle  which  is  not  sharp  edged,  such  that  the  separation 
line  moves  with  change  in  velocity,  the  Reynolds  number  affects  the  drag 
coefficient,  and  compensations  will  have  to  be  made  for  this  effect.  A 
po*sibility  exists  in  artificially  tripping  the  boundary  layer  on  the 


obstruction  so  as  to  induce  turbulence  locally  and  fix  the  boundary 
layer  separation  line.  However,  such  refinements  have  not  been  used  in 
this  study,  which  is  intended  to  furnish  qualitative  information  rather 
than  quantitative  design  data  and,  in  that  case,  it  is  unnecessary  to 
substantiate  the  small  improvements  in  similarity  which  can  be  had  by 
artificially  inducing  separation  on  the  model  hill.  Thus,  the  problem 
of  scaling  need  not  concern  us  in  this  study,  especially  since  a  com¬ 
parison  with  field  data  is  not  possible  at  this  time.  We  shall,  there¬ 
fore,  formulate  our  problem  in  more  detail  without  regard  to  scaling. 


Chapter  il 

THEORETICAL  CONSIDERATIONS 


The  two  essentially  different  problems  which  arise  m  considering 
the  interaction  of  missiles  and  turbulent  wind  fields  are  that  of  missile 
flight  stability,  and  that  of  impact  dispersion.  The  difference  of  these 
two  problems  can  best  be  illustrated  by  considering  the  flight  of  a  missile 
through  a  homogeneous  velocity  field  of  infinite  extent.  A  missile  which 
flies  at  constant  speed  encounters  a  spatially  random  velocity  field  which 
is,  with  respect  to  a  coordinate  system  traveling  with  the  missile,  con¬ 
verted  into  a  random  and  stationary  time  series  of  the  continuous  variable: 
velocity.  If  the  missile  has  a  transfer  function  H(w)  ,  then  the  missile 
response  velocity  spectrum  <p  (w)  is  related  to  the  impact  wind-gust 
spectrum  $  (u>)  by  the  relation  of 

<t>  _,(<*>)  = 
m 

Thus,  since  the  transfer  f'-ncuion  {  H(u>)  f  is  a  deterministic  function, 
and  since  <j>  (co)  for  an  infinitely  long  stationary  record  denotes  the 
exact  average  behavior  of  the  wind  field,  <p  (<■>)  is  also  an  exact  average 
measure  of  the  missile  response.  If  none  ofthe  response  amplitudes  ex¬ 
ceed  the  stability  limit  of  the  missile,  then  only  some  fluctuations  of 
the  missile  occur;  if  some  do  exceed  the  stability  limit,  the  missile 
might  change  course  drastically  and  miss  its  target  by  a  wide  margin. 

The  s. ability  can  usually  be  evaluated  on  the  basis  of  the  average  be¬ 
havior  expressed  by  Eq.  2.  In  this  paper,  we  shall  provide  experimental 
data  on  wind  spectra,  which  can  be  used  for  missile  stability  calculation 
purposes . 

In  contrast  to  stability,  the  dispersion  of  a  missile  results  from 
an  integrated  effect  of  all  the  velocities  which  are  acting  on  the  missile 
in  its  course  along  the  missile  trajectory.  Since  these  velocities  are 
fluctuating  from  instant  to  instant,  and  can  be  described  only  in  a  proba¬ 
bilistic  way,  the  missile  dispersion  cannot  be  predicted  deterministically. 
Instead,  the  missile  dispersion  problem  is  the  problem  of  determining  the 
probability  distribution  of  the  missile  trajectory  end  point  as  a  function 
of  the  sequence  of  all  the  velocities  which  the  missile  encountered  along 
the  trajectory.  The  distribution  of  the  end  point  of  the  missile  then 
becomes  a  function  of  the  joint  probability  distribution  for  all  the 
velocities  along  the  missile  trajectory. 

In  this  report,  we  shall  disregard  the  characteristics  of  the 
missile  and  shall  concentrate  on  an  attempt  to  describe  the  joint  proba¬ 
bility  distribution  for  the  velocities  along  some  hypothet ical  missile 
trajectories  in  a  simplified  manner.  The  theoretical  ideas  will  be 
developed  in  this  chapter.  They  lead  to  a  p  jgram  of  measurements  of 
probability  distributions  which  was  performed  in  the  Fluid  Mechanics 
Laboratory  of  Colorado  State  University. 


H(u>) 


(<*>) 

w 


(d) 
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Since  it  is  impossible  to  obtain  the  true  joint  probability  dis¬ 
tribution  for  all  velocity  vectors  along  any  trajectory,  a  simplifying 
procedure  has  to  be  adopted.  We  proceed  by  introducing  some  simplifying 
assumptions  wnui  represent  the  turbulence  encountered  by  the  missile 
by  the  instantaneous  turbulence  existing  along  the  mean  trajectory. 

Furthermore,  the  trajectory  is  subdivided  into  sections  and  it  is  assumed 
that  the  turbulence  in  each  section  can  be  represented  by  the  turbulence 
at  the  end  points  of  the  sections.  For  the  ensuing  sequence  of  velocities 
at  the  section  end  joints,  the  joint  probability  density  function  is  then 
cons  >.  :uc red  and  broken  down  into  a  product  of  functions  which  can  be 
determined  by  means  of  available  experimental  techniques.  No  attempt 
will  be  made  to  apply  the  ensuing  functions  to  the  missile  dispersion 
problem. 

1 

2 . 1  Bas i c  Assumptions  j 

The  problem  of  evaluating  the  instantaneous  missile  trajectory  is  S 

approached  in  the  following  way.  Let  the  mean  trajectory  of  a  missile  be  ) 

given,  and  use  the  reference  coordinate  system  as  shown  in  Fig.  1  ft"*  our  i 

problem.  Then  on  its  travel  along  the  trajectory  the  missile  encounters 
mean  velocities  and  a  sequence  of  gusts,  both  described^by  a  velocity  \ 

vector  v(s;t)  ,  where  t  is  the  time  of  flight,  and  s  is  the  position  j 

vector  of  the  trajectory.  The  velocity  Rector  consists  of  a  mean  vlocity 
V(s)  and  a  fluctuation  in  velocity  v'(s;t)  .  The  position  vector  con¬ 

sists  of  a  mean  position  vector  J  corresponding  to  an  absence  of  all 
velocity  fluctuations  U-e-,  the  trajectory'  due  to  mean  wind  only!  and 
a  small  deviation  ■;  -s  due  to  the  sequence  of  fluctuating  velocities 
which  the  missile  has  encountered  during  the  time  t 

Now,  let  the  travel  time  until  impact  be  equal  to  t.  and  the  end 
point  of  the  mean  wind  trajectory  be  located  at  x  .  T^ien  due  to  the 
sequence  of  wind  fluctuations  encountered  during  its  flight,  the  missile 
is  deflected  in  the  impact  area  by  a  total  deviation  r'  from  the  target 
distance  x  Due  to  the  random  nature  of  the  fluctuations  encountered, 

the  r'  will  also  be  randoml)  distributed.  The  probability  distribution 
of  the  quantity  r'  is  the  desired  quantity  to  which  the  results  of  this 
study  must  be  applied. 

The  meteorological  problem  associated  with  finding  the  probability 
distribution  of  v\  is  to  make  available  knowledge  of  ■.!»«  instantaneous 
velocity  field  which  the  missile  might  encounter  on  its  course.  Clearly, 
this  problem  cannot  be  solved  by  presently  available  techniques.  In¬ 
stead,  it  is  proposed  to  obtain  joint  probability  distributions  for  the 
simultaneous  occurrence  of  a  sequence  of  velocity  vectors  along  the 
missile  trajectory.  In  general,  this  requires  specifying  joint  proba¬ 
bility  distributions  of  the  ’oint  occurrence  of  velocities  at  infinitely 
many  different  points  in  space  and  time.  In  order  to  reduce  this  problem 
to  tractable  dimensions,  a  number  of  assumptions  have  to  be  made. 

The  first  assumption  is  that  the  distance  of  any  instantaneous 
trajectory  from  the  mean  trajectory  calculated  on  the  basis  of  the 
mean  wind  distribution  is  small,  so  that 


v '  ( s  ;  t  )  ~  v  ’  (  s  ;  t ) 


() 


In  this  manner,  it  is  no  longer  necessary  to  consider  the  whole  space 
but  one  can  concentrate  on  the  single  trajectory.  Obviously,  the  validit 
of  this  assumption  depends  both  on  the  relative  magnitude  of  v'  with 
respect  to  v  ,  and  on  the  response  characteristics  of  the  missile,  and 
will  have  to  be  tested  each  time. 

The  second  assumption  concerns  the  time  distribution,  he  assume 
that  the  missile  travels  much  faster  than  the  velocity  fluctuates,  so 
that 


V  (s  ;t)  v ' (s  ;t  ) 


where  t  denotes  the  start  time.  This  assumption  implies  that  during 
the  flipRt  time  the  relation  holds: 


v ' (s  ;  t )  v'(s;to)  pj  v'(s;tQ)  v '  ( s  ;  t  q  ) 
or  that,  in  the  average  for  n  different  starting  times  t 


n 


1  . 


n 


E  v'  S;t  v'(3't  .  -  -  ::  (,v  (S ; t  .  )  )■ 

.  ,  ’  01  n  .  ,  01 

i=l  i=l 


If  the  flow  is  stationary,  and  if  the  ergouic  hypothesis  is  valid,  then 
we  can  restate  this  requirement  as: 


R  ~  1 

T 


where 


R  is  the  autocorrelation  function  defined  bv : 


,  T  v  '  (  P  ;  t  *  !  t  - 1  1  )  •  v  '  1 5  ;  t  ^  d  t 
R  -  R(t  -t  .  =  i  f - ^ - - 2. 

T  X  0  id 

o  ^ 

V '  ■  [s  ,  t  1 
0 


'  1 1 


whe  re 


T  is  an  observation  time  t  awen  long  enough  to  ensure  a  stable 
average,  and 

t  -t  is  the  1 1  me  dti-ing  which  the  missile  has  traveled  from 

0  X  t  O  X  * 

o 


To. convert  actual „t  ravel  times  to  model  travel  t  ires  t'.e  »cil 
.  t  u  * ,  f  u  . 

. — — 1  =  ! - 1  must  fe  used,  wh l eh ,  tor  u*  : 


“  o  mode  1 
reduces  to 


~ o  field 
' field 


mode  1 


.  i .  s  . .  i .. 

f  ;  e  i d 


o  mode  1 


"o  f’.eh 
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For  small  times  t  -t  Eq.  4a  becomes: 
x  o 


(t  -t  ) 
X  o 


(4b) 


V 


whe  re 


is  the  microscale  of  the  turbulence.  The  scale  X  can  be 


replaced  to  a  good  approximation  by  the  scale  X  of  the  u-component 
of  the  turbulence 


X'2 

g 


(5) 


Consequently,  it  follows  that 
to  be  valid. 


tx-tQ  <<  Xg  for  the  assumption  Eq.  3 


2 . 2  Simplifications  of  the  Probabilistic  Problem:  connecting  proba¬ 
bilities  along  trajectories.  ” 

We  base  our  calculations  on  assuiqptions  Eq.  2  and  Eq.  3,  and,  thus, 
we  have  reduced  the  meteorological  aspects  of  the  problem  to  finding 
simultaneous  instantaneous  velocity  distributions  along  the  mean  tra¬ 
jectory  x  .  To  avoid  the  implied  necessity  of  determining  velocities 
simultaneously  at  infinitely  many  different  points,  we  adopt  the  follow¬ 
ing  probabilistic  specification  of  the  velocity  field.  The  required 
quantity  is  the  joint  probability  density  distribution 


f. 

J 


(6) 


for  all  n  points  along  the  mean  trajectory.  The  experimental  distri¬ 
bution  of  f  requires  simultaneous  measurements  at  all  n  points  of 

the  trajectory,  i.e.,  it  requires  an  infinite  array  of  probes  placed  along 
the  trajectory.  Evidently,  this  is  an  impossible  task,  so  that  instead, 
the  trajectory  is  cut  into  n  finite  intervals,  of  length  Ax  ,  at 
whose  end  points  turbulent  quantities  are  measured.  In  each  interval 
Ax  =  xi+i-xi  the  instantaneous  velocity  is  assumed  constant  and  equal 


v.'  =  u.'i  +  v.'  i  +  w ! ic  (7) 

i  i  i-  i  ’ 

when  the  components  u.'  ,  v.'  and  w?  are  average  values  of  the  velocity 

components  at  the  two  Ind  points .  From  the  values  of  v.'  ,  the  tra¬ 

jectory  is  calculated. 

The  problem  to  be  solved  then  is  to  convert  probability  distri¬ 
butions  between  adjacent  points  in  such  a  way  that  a  meaningful  approxi¬ 
mation  for  Eq.  6  is  found.  We  want  to  investigate  three  smiple  cases 
of  possible  approximations  for  Eq.  6. 


8 


a.  Consider  first  the  assumption  that  f(v*  )  and  f(vJ) 
are  statistically  independent.  This  condition  corresponds  to  velocities 
which  vary  comparatively  rapidly  along  the  trajectory,  in  the  sense  that 
R  ^  0  where  R  is  the  spatial  correlation  coefficient  obtained  from 

v  X  r 

tne  definition 

R  .  R„.  ,-*.)  .  i-  f f 

x  1+1  1  Ax  J  , . — - - -  =====  .  (8) 

Xi  */v,2(tQjxi)  •  /v,?  (tQ,xi  +  1) 

However,  the  assumption  of  rapidly  varying  velocities  is  in  contradic¬ 
tion  to  the  assumption  of  a  velocity  vector  which  is  constant  throughout 
the  travel  interval  Ax  ,  unless  Ax  is  chosen  in  such  a  way  that  a 
meaningful  relation  between  it  and  the  space  integral  scale  *Y  exists, 
where : 


•  J.  R, 


Also,  in  order  to  be  of  influence  on  the  flight  pattern,  'J  must  be 
large  compared  to  the  length  dimension  L  of  the  missile,  such  that  a 
condition  for  the  validity  of  this  assumption  might  be  defined  as: 


ax~  ys 


<<  1  say  <0.1 


Under  these  circumstances,  Eq.  6  reduces  to 


f(v'Q  ,  vj  ....  W  j  ,  vp  =  f(*»)  f(vj) . fCv^j)  f(v^)  (11a) 


or  in  terms  of  conditional  probability  densities: 
fCv;  I  ,  v;_2 . )=  tlfp 


(lib) 


This  equation  can  be  evaluated  conveniently,  if  the  probability  density 
distributions  f(v<)  are  given.  These  correspond  to  joint  probability 
densities  for  thos£  variables  uf  ,  vf  and  w/  ,  which  will  be 
discussed  below  in  Section  2.3.  1  1  1 

b.  As  a  second  possibility,  we  considered  the  condition 


Js  > 


x  -x 
e  o 


in  which  case  the  correlation  coefficient  defined  by  Eq.  8  assumes  a 
value  very  near  to  1.  This  implies  that  the  velocities  v'(t  ,  x.) 
v'(t  ,  x.  +j)  are  very  nearly  proportional,  so  that  0  1 


v*  (t  ,  x.  , )  js  Av'  (t  ,  x. ) 
o’  i+l  o’  1 


9 


where  A  is  a  (vector)  constant.  Furthermore,  the  jpdf  defined  by 
Eq.  6  becomes: 


f(v;) 


or  in  terms  of  conditional  probability  densitir 


(13a) 


f(v'  i  v!_!  ,  v?_2  . )  =  1  (13b) 

Again,  the  discussion  of  a  method  for  calculating  f(v^)  is  postponed 
until  Section  2.3. 

c.  The  assumption  of  a  and  b  bracket  the  possibilities  for 
..simplifying  the  joint  probability  density  functions  of  the  turbulence 
along  the  trajectory.  An  intermediate  method,  based  on  the  assumption 
that  the  eddy  structure  of  the  turbulence  is  highly  elongated,  (as  is 
usually  the  case  in  turbulent  flows)  would  combine  assumptions  of  in¬ 
dependence  of  the  motions  perpendicular  to  the  mean  wind  direction  with 
an  assumption  of  some  dependency  of  the  components  in  the  wind  direction 
along  the  trajectory.  The  simplest  way  is  obtained  if  a  Markoff  de¬ 
pendency  can  be  found  to  relate  probability  density  distributions  along 
the  trajectory,  i.e.,  if 

f(v!  |  v^_x  ,  v;_2 . v{,  vp  =  f(v*  |  *U)  (14) 

when  f (v!  |  v!  ^)  is  the  conditional  probability  density  for  the 

occurrence  of  v!  when  v!  ,  has  already  occurred. 
i  l-l  J 

The  elongated  eddy  structure  leads  us  to  assume  that  what  happens 

at  point  depends  on  the  happenings  at  ^  only  through  the  ul- 

component ,  i.e.,  the  components  v!  and  w!  are  independent  of  all 

components  at  the  point  x^  ^  except  inasmuch  as  they  depend  on  u!  , 

which  in  turn  depends  only  on  the  component  u!  .  and  not  on  v.  .  and 

w.  .  .  Write  F.q.  14  in  the  form:  '  1 

l-l  n 

f(v|  [  v!_j)  =  f(u!  |  u|_x)  f (v!  |  u!)  •  f(wl  |  u!,  v!)  (15) 


where  f(w.'  |  u!,  v? )  denotes  the  conditional  probability  density  for 

finding  wf  when  both  u!  and  v!  are  assumed  to  occur  also. 

*  i  ii 

We  can  now  summarize  the  results  for  the  three  approximations  of 
Eq.  6  as  follows: 


Independence  (Case  a) 


) 


so  that 

ffV  vi-l’  Vi-2 ' 


,v;)  = 


f(v!) 


f(v?_ 


•f(V;) 
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Dependence  (Case  b) 


v 


i-1  * 


)  =  1 


so  that 


f(v. 


,)  = 


f(v.) 


Markoff  dependency  (Case  c) 


f  ( v  .* 


C-r  v i - 2 ’ . 1  s  f(vi  l  vl-i) 


so  that 

f(v.,  v.  ,,  v.  . .v  )  =  f(v’)  f(v'  !  v'J . f(v?  1  v?  ,) 

1  *  1-1  ’  1-2  o  0  1  0  1  1-1 

whicn  simplifies  further  for  the  elongated  eddy  case  to  Eq,  15. 


2 . 3  Simplification  of  the  Probabilistic  Problem:  joint  probability 
densities  at  a  point 


All  three  cases  discussed  above  require  the  determination  <jf 
probability  density  distributions  of  the  form:  ^f(v?)  .  Since  v! 
is  a  vector  consisting  of  three  components,  f(vp  1S  actually  a 
joint  probability  density  function  for  the  joint  occurrence  of  uj  , 
v!  and  .v!  .  Such  a  triple  joint  probability  density  function  is  too 

difficult  io  determine  experimentally.  We,  therefore,  write 


f  (u! ,  v.' ,  w! )  =  f  (u! ) 

v  l  ’  i  ’  i  l 

in  the  form 

f (u!  v.',  w!)  =  f(u') 

1*  1  1  1 

and  Eq.  14  in  the  form 

f(vj  |  v ! )  =  f (u!  | 


f(v!  |  u!)  f(w^  j  u!,  v?) 
f(v|  I  u!)  f(w') 
u!  _1]  f(v!  |  up  f(w!) 


(16) 


(17) 


(IS) 


which  are  based  on  the  following  assumptions: 

a.  The  velocity  component  w.'  is  statistically  independent  oc 
all  other  velocity  components. 

b.  The  connection  between  adjacent  points  takes  place  only  through 

u!  and  is  at  most  first  order  Markovian. 
i 

Assumption  a.  is  partly  justified  because  the  homogeneity  of  the 
turbulence  in  planes  parallel  to  the  ground,  in  a  two-dimensional  flow 
field,  requires  that  the  time  average  product  u|  wj  =  0  which  is  a 
necessarv,  but  not  a  sufficient  condition  for  statistical  independence. 


1 1 


The  fi rsr  part  of  assumption  b.  is  postulated  without  any  firm 
basis  except  for  the  motion  of  an  elongated  eddy  stated  previously. 

For  the  second  part,  however,  we  have  some  support,  both  from  meteor¬ 
ological  data  as  well  as  for  the  laboratory  case  of  the  present  study. 

For  a  Markoff  dependency  to  exist,  it  is  a  necessary  and  sufficient 
condition  that,  if  the  variables  u{  are  stationary  with  respect  tc  i 
and  Gaussian,  and  are  a’so  jointly  6aus.si.an  distributed,  then  the  cross  j 

correlation  between  u.'  ,  u.'  +  J  ,  u'  ^  is  an  exponential  function  j 

in  the  parameter  i  (p.  96,  Feller  (1964)  p.  234,  Doob  (13531 li.e.,  in  j 

the  continuous  parameter  case  5 

R(x)  -  e  vXR(Oj  x  0  (19)  | 

s 

where  A  has  a  non -negative  real  part,  if  R(x)  is  known  to  be  con-  | 

tinuous.  Conversely,  a  sequence  with  stationary  Gaussian  distributions  5 

satisfying  Cq.  18  is  Markovian  and  Hq .  14  can  thus  be  used.  f 

t 

i 

The  applicability  of  a  Markoff  process  to  turbulence  data  is  thus 
insured  if  it  can  be  shown  that  i 

a.  the  space  correlations  are  homogeneous,  i.t-.,  independent 

from  where  the  correlation  starts.  • 

b.  the  space  correlations  are  exponential. 

c.  the  probability  density  distributions  for  the  functions  repre¬ 
senting  u'  are  Gaussian. 

Some  proof  for  the  validity  of  these  conditions  for  our  laboratory  flows 
will  be  given  in  the  next  chapter,  for  atmospheric  turbulence  in  neu¬ 
trally  stratified  atmospheric  boundary  layers  over  homogeneous  terrain, 
these  conditions  are  approximately  satisfied.  In  the  older  meteorological 
literature  (for  reference  see  Pasauill  (1961))  the  autocorrelation  functions 
were  usual’.-  found  to  be  exponential.  Together  with  Taylor's  hypothesis, 
according  to  which  time  correlations  can  be  translated  into  space  corre¬ 
lations  by  means  of  the  substitution  t  =  —  ,  (where  0  is  the  local 

mean  velocity),  it  can  thus  be  shown  that  space  correlations  are  expo¬ 
nential.  Meteorologists  have  in  recent  times  (Lumley  and  Panofsky  (1964)) 
preferred  to  use  different  analytical  representations  of  the  correlation 
functions,  for  the  simple  reason  that  the  spectrum  corresponding  to  an 
exponential  autocorrelation  decreases  at  large  values  of  angular  fre¬ 
quencies  w  proportional  to  uT;  ,  whereas  the  spectral  shape  should 
contain  an  inertial  subrange.  With  a  drop-off  proportional  to 
The  difference  between  5/3  and  2  is,  however,  not  large  enough  to 
give  a  strong  reason  for  discarding  the  assumption  of  an  exponential 
decay  of  the  autocorrelation  function.  For  our  prediction  purposes,  it 
docs,  therefore,  seem  to  be  justified  to  assume  that  an  exponential  auto¬ 
correlation  fun Cl  ion  exists  in  neutrally  stratified  atmospheric  boundary 
layers  over  homogeneous  terrain.  In  a  later  section  we  shall  show  that 
an  approximately  exponential  space  correlation  function  which  is  homo¬ 
geneous  along  trajectories  parallel  to  the  ground  can  be  found  even  in 
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the  highly  disturbed  flow  field  downwind  from  a  mode)  ridge.  Since  wc 
also  find  that  almost  all  velocity  components  follow  Gaussian  distribu¬ 
tions,  the  Markoff  dependency  postulated  for  tq .  15  is  reasonably  well 
established  experimentally. 

2.4  Some  Considerations  on  Gaussian  Two  Variable  Joint  Probability 
Dens  i  tv  Tunct  ions  „ 

— - i. — _ ___ 


When  the  joint  probability  density  functions  of  the  quantities  of 
turbulence  at  one  point  are  Gaussian,  then  this  distribution  function  is 
fully  specified^ by  the  means  and  the  turbulence  quantities  u'2  =  a  -  , 

v 1  ^  =  c  •’  ,  w'-  =  •  •  as  well  as  by  the  cross  correlations,  for  examp] 

u’v'  .  These  quantities  are  most  important  also  in  describing  the  dy¬ 
namic  conditions  of  the  turbulent  flow,  i.e.,  they  represent  stresses, 
and  it  is,  therefore,  of  interest  to  show  the  connection  between  the 
probability  density  functions  and  the  stresses. 


Theoretically,  if  all  the  probability  densities  of  individual 
turbulent  components  are  distributed  in  a  Gaussian  fuii.i,  then: 


f(u')  = 


>2-  c 


u' 


■  (u' -m  , V 
_ u' 

2  0 

u' 


f(v')  = 


/2n  c 


2a  ,2 

v' 


(20) 


-(w'-m,,)' 


f(w')  = 


t/Tn 


2a 


The  joint  densities  of  two  turbuient  components  can  be  exnressed  by  a 
joint  Gaussian  form,  i.e., 


f(u'.v')  * 


-Q(u',v') 


2"  a 


/l-c 


(21) 


for  some  constants  a  ,  •  0  .  a  ,  >  0 

u'  v'  ’ 

<  m  ,  <  +<®  ,  in  which  the  function  Q(u', 

number^  u?  and  v'  is  defined  by 


o  < 

V’) 


for 


<  m  ,  <  +® 

u' 

any  two  real 


13 


1  =  1 

V' 

u 1  - m  .) 

j  v'-m  , 

/ 

v'-m 

->  i  u'  1 

1  V 

l 

♦  1 

V  '  1  1 

.  i  i 

2(l-p?) 

LI"., 

\  u'  / 

"  -1  f 

\  V  / 

i  °v 

l 

where  r  is  the  correlation  coefficient 


u'v ' 


su7?  fy7?- 


and 


,  are  the  variances  of  u'  and  v 1  .  resDective  1  v .  and  m  ,  and 
m_,  are  the  mean  values.  The  curve  Qiu'.,  v';  =  constant  is  an  ellipse 
since  p  1 

In  order  to  find  the  orientation  of  the  ellipse,  the  coordinates 
u*  and  v*  of  the  coordinate  system  parallel  to  the  axes  of  the 
el  1  ipse  • 


u*  =  u'  cos  *'  +  v1  sm 


(22) 


v*  =  -u'  sin  •;  +  v’  cos 


5  r»  o  q  A 


are  introduced.  Applying  the  Jacobian  transformation  to  the  probability 
density,  Fq.  20,  we  obtain 


f(u*,  v*)  = 


f  ( u ' ,  v ' ) 


!j(u',v')| 

=  f(u*cos  6  -  v*  sin  6,  u*  sin  0  +  v*  cos  8) 
which  follows  from  the  fact  that  the  Jacobian: 


(23) 


j(u',v')  = 


au* 

au* 

au' 

a  v ' 

ov* 

ov* 

Q  1 

C 

ov' 

=  1 


For  further  simplicity,  we  mav  assume  m  ,=  m  ,=  0  then, 
‘  •  u’  v' 


f(u*,v*)  = 


2"o  ,  o  /1-p" 

u '  V  ' 


exp  <  - 


1 


2(l-p2) 


-2p 


\  u' 


0  . 

u'  v' 


s  ur 


o 

v' 


U*' 


(24) 


cost1  Sinr  1  S II. 


cost)  sine 


u*v ' 


u' 


0,0. 

u'  v' 


V  1  \ 


/ 


|  sin-'e  cost)  sinf  cos-’f 

+ - +  2p - +  -  ! 


,  I 


o 


U  f  V 


i 
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because  of  the  symmetry  of  the  ellipse  with  respect  to  the  new  axes,  the 
term  involving  u*  v*  should  vanish,  i.e., 


cos8  sine 


-P 


sin20-  cos2e 


cos0  sine 


=  0 


This  can  be  written  in  terms  of  the  turbulent  stresses  by  means  of  the 
relations : 


and 


o 


u' 


tan  26 


2u-v' _ 

u*2-  v'2 


(26) 


With  this  relationship,  the  joint  density  function  of  the  turbulent  com- 
pc-’ents  can  be  defined  once  we  have  the  values  of  the  associated  turbu¬ 
lent  stresses.  For  example,  the  Eq,  18  can  be  established  by  measuring 

u!2,  ,  u!2  ,  v!2  ,  w.'2  ,  u?  v?  and  u!  ,u!  .  For  its  application, 

we  have,  however,  to  shew  that  the  individual  components  are  Gaussianly 
distributed  and  that  the  joint  probability  distributions  follow  Eq.  19. 


We  notice  in  passing  the  equality  between  Eq.  26  and  the  inclination 
of  the  plane  of_&ero  shear  stress  in  a  plane  stress  state  of  classical 
mechanics,  if  u'2  and  v*2  are  the  normal  stresses  and  u'v'  is  the 
shearing  stress  (due  to  turbulence).  Clearly,  then,  the  angle  9  denotes 
the  orientation  of  a  plane  at  a  point  in  a  fluid  when  the  shear  stress  is 
zero,  so  that  the  joint  probability  distribution  is  found  to  be  oriented 
with  the  long  axis  of  the  ellipse  of  constant  correlation  parallel  to  the 
zero  shear  stress  plane. 
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Chapter  III 

EQUIPMENT  AND  PROCEDURES 


^ . 1  The  Experimental  Setup 


The  experiments  were  performed  in  the  U.S.  Army  Meteorological 
Kind  Tunnel  in  the  Fluid  Dynamics  and  Diffusion  Laboratory  of  Colorado 
State  University.  This  facility  is  shown  in  Fig.  2.  It  is  a  recir¬ 
culating  wind  tunnel  with  an  88  ft  long  test  section  with  a  6  x  6  ft2 
cross  section.  For  the  experiments  of  this  study,  the  model  hill  was 
placed  at  a  distance  of  approximately  40  ft  downstream  from  the  inlet 
where  the  undisturbed  boundary  layer,  stimulated  by  large  roughness 
elements  in  the  inlet  region  of  the  test  section,  had  an  undisturbed 
thickness  of  about  24  inches.  The  model  hill  consisted  of  a  plexiglass 
section  with  a  shape  n  given  by 


v 


h 


x 


COS  -j — 


for 


i 

•> 


x  <  1^ 

L  "  2 


(27) 


where  the  base  width  L  =  20  in.  and  the  height  h  =  4  in.  Th"  velocity 
outside  of  the  undisturbed  boundary  layer  was  30  fps . 

3.2  Measurement  of  Mean  Velocity  Profiles 


Mean  velocity  profiles  were  measured  both  by  hot  wire  anemometer 
and  pitot  tube,  in  order  to  obtain  a  cross  check.  In  the  upper  part  of 
the  flow,  continuous  traverses  of  velocity  were  taken.  In  the  lower  part 
or  in  the  separation  region  where  the  variability  of  velocity  was  large, 
point  by  point  data  were  taken  in  order  to  determine  the  velocity  pro¬ 
files  more  precisely. 

The  hot  wire  measurement  of  mean  velocity  was  made  with  a  4  x  10  b 
inch  diameter  single  wire  which  was  held  perpendicular  to  the  local  mean 
velocity  vector  0  .  The  hot  wire  anemometer  used  was  made  at  CSU. 

By  means  of  the  pitot  tube,  total  head  readings  were  obtained  for 
calculating  mean  velocities.  If  there  is  no  pressure  gradient  in  the 
flow  field,  the  jeal  mean  velocity  can  be  calculated  by 


1 

1 


pO2 

a 


AP 


AB 


(28) 


AP^g  =  pressure  difference  between  the  static  tap  and 
dynamic  tap  of  a  pitot  static  tube. 

But  in  th neighborhood  of  the  model  ridge,  large  pressure  gradients 
exist,  not  only  in  y-direction  but  also  in  x-direction.  Since  the 
static  tap  is  one  inch  downstream  from  the  dynamic  tap  on  the  pitot 
tube,  a  correction  must  be  applied  for  the  pressure  gradient  between 
the  twro  taps.  Since 


16 


7  ?2  =  iPAB  -  <V  V  <-’9) 

p  =*  density  of  air  at  the  room  temperature 

a 

U  =  local  mean  velocity 
APab  =  the  measured  pressure  difference 

-  the  static  pressure  at  dynamic  tap's  position 

Pg  =  the  static  pressure  at  static  tap's  position 

P  -  Pg  =  pressure  difference  between  the  static  tap  and  the 
dynamic  tap. 

If  p  -  p  is  know.i,  the  local  mean  velocity  at  one  point  can  be  cal- 

A  B 

culated  from  Lq.  29.  At  each  point  the  value  of  PA  -  Pg  can  be  ob¬ 
tained  from  Fig.  3.  This  figure  was  made  by  connecting  the  static  tap 
and  a  reference  tap  to  the  pressure  transducer  (Transonic  Type  120  Equibar) . 
Since  the  static  tap  is  one  inch  downstream  from  the  dynamic  tap,  at  one 
point  the  coordinates  of  the  dynamic  tap  is  known,  say  (x.,  y^)  then 
the  static  tap  is  (Xj  +  1,  y^  .  When  the  coordinates  of  two  points  are 

known  the  pressure  difference  P  -  Pg  can  be  obtained  from  Fig.  3  and 

the  corrected  mean  velocity  at  that  point  can  be  calculated  by  applying 
Eq.  29. 

To  measure  the  mean  velocity  profiles  the  pitot  tube  and  the  hot 
wire  were  mounted  on  a  24  inches  vertical  carriage.  The  dynamic  tap  of 
pitot  tube  and  hot  wire  were  held  side  by  side  at  the  s«me  height.  The 
velocity  profiles  were  taken  every  two  inches  downstream  from  the  crest 
up  to  x  *  18"  and  also  at  x  =  24",  36",  40" 

The  block  diagram  of  set  up  is  shown  on  Fig.  4.  The  calculation  of 
tiie  mean  velocity  is  listed  on  Table  I  and  the  results  are  on  Fig.  5. 

When  a  hot  wire  was  used  to  measure  the  mean  velocity,  the  cali¬ 
bration  curve  of  this  wire  was  checked  from  time  to  time  and  the  wire  was 
recalibrated  if  excessive  drift  of  an  anemometer  was  detected.  It  was 
found  that  after  a  hot  wire  had  aged  several  hours  the  drift  of  the  wire 
was  negligibly  small. 

3.3  Measurement  of  Turbulent  Quantities 

For  coordinates  of  the  flow  field  as  shown  in  Fig.  1,  the  turbulent 
components  at  a  point  in  x,  y,  z  directions  are  u',  v',  w  * ,  respectively. 

The  u' -component  was  calculated  from  a  single  hot  wire  held  parallel 
to  2  axis  (Fig,  61.  The  v' -component  was  calculated  from  a  crossed 
wire  held  in  the  x-y  plane  (Fig.  6). 
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In  the  subsequent  discussion,  we  shall  use  the  following  notation: 


u  = 


v 


and 


w  =  /w'‘-  is  the  rms  value  of  the  fluctuating  velocity  component 
in  the  x,  y,  z  direction,  respectively 

uv  =  covariance  of  the  fluctuations  u'  and  v' 


e 


1, 


are  the  rms  values  of  ^he  fluctuating 
voltages  ej  and  e^  measured  with 

wire  No,  1,  or  2,  respectively. 


a.  Calculating  of  u-component  of  turbulence. 

To  calculate  the  u-component  at  one  point,  savr,  (x^,  y^)  we  need 
the  following  information: 

1.  the  rms  value  of  a  single  wire  at  (x,  y) , 

2.  the  calibration  curve  of  this  wire, 

3.  the  local  mean  velocity  U  at  (x,  y) 

dE  -r 

4.  the  slope  of  the  calibration  curve  at  0  then: 

dE 

°  3U  U 


I 

1 


) 


or 


I  dll  I  •’  • 

u‘  =  m  e‘ 

The  u-component  at  a  point  (x,  yl  was  calibrated  by  Eq.  30. 
b.  Calculating  v  and  uv  components  of  turbulence. 


(30) 


If  the  crossed  wire  is  held  in  the  x-y  plane  as  shown  in  Fig.  7, 
then  we  find  in  general  that 

1  - — — 

ci  =  -7ZTT-  dir  ,(u’cos  1  +  v'sln  •■'')- 

(31) 


e ,  -  — — ; — —  2  >  i u ’ c os  ;•  -  v’sin  a )  • 

k  C  O  5  i*  fr-; - 

JU 

dlij  dl  : 

where  ^p  and  are  the  slopes  of  the  cal  i  brat  ion  urve  for  wire  l 

and  w  i  rc  .  respective!)  ,  Equation  31  can  be  written  as 


I  dE  \  1 


1  ■> 

COS"  a 


1  ~~o 

(ti'^cos  a  +  C'".,v'  cos  a  sin  a  +  v'*  sin’al 


.  dE  \2  - 

?  1  |  l  l  2  - *  - - 

e~  = - i  -T77-  (u'‘cos  c  +  au'v'  cos  8  sin  8  +  v'-  sin-2) 

2  du 

cos“B  \ 

In  order  to  account  for  small  deviations  of  mean  velocity  vectors  from 

the  horizontal  we  write: 

,,o 

a  =  45  +  i 

8  =  45°  -  i 

where  41  is  the  small  deviation  of  the  angle  between  the  mean  velocity 
vector  and  the  horizontal,  as  determined  from  the  streamline  pattern . 
Then,  for  small  4  ,  so  that  rns  in  can  be  neglected,  and 

cos 4  ~  1  ,  s  iny s  ;  : 


cos  a  *  cos  (45  ♦  t )  =  -r- 


11-;) 


cos^n  = 


(!♦:) 


s  in*  a  =  *•  ♦ 


cos  r  =  (  1  *: ) 

i 

cos*  c  =  ^  +  :■ 

.  -n 

►  C 

1 

s  1  n  h  =  T  ( 1  - ;  i 

s  1  n  ‘  r  -  —  -  t 

Substituting  !:q  ,  55  1 

nto  i.u.  52  yields 

/  ' 

.JdU 

l(y  •  i 


e 2 ( ■y  '  ;)i  tt"  “  ( t  -  ; )  *  u 1  v '  +  v 1  ■  i \ 


1  ,  —  1 

u-  (T  -  ;  1  *  uv  ♦  r 


i +  iarj 


; !  ♦  u’,‘  -  v*  (4  -i-l 


The  difference  of  Eq ,  54  and  !.q. 


1  I  dll 


.v>  tnen  vieids 


1  dl) 


■  iu-  -  v 
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Inserting  Eq.  36  into  Eq.  34  leads  to 


2  1 
v  =  7 


/ 


dU 

dE, 


/ 


+  e. 


dU 

3E7 1 


-  u 


(37) 


r  ' 
«  2 

1 

dU 

2  „  2 

-i  - 

dU 

2 

61 

\^i 

*  e2 

dE2  J 

Since,  in  our  experiment  both  wires  were  of  the  same  length,  the  slope  of 
the  calibration  curves  of  both  wires  for  the  same  U  wire  the  same,  i.e., 

at  a  point  .  Then  Eq.  37  becomes 


2  1  dU 

1  ■  7  JE 


and  Eq.  33  yields : 


(e^2  +  e^2)  -  u2  -  4> 


dU 


'<• 1 2  -  e2  ) 


(38) 


u'v'  =  J  ^  Cj2  Cj  -  4>)  -  e22  (j  +  $)  +  <J> (u2  -  v2)  (39) 

Substituting  hq.  38  into  Eq.  39 


—  1  ,  2  (JU?  ,  ,  ? 

uv  »  J  (Cj‘-  e22)  ^  -  2  *v2 


dE  "  ‘  T’  (40) 

If  <j>  =  0  i.e.,  the  velocity  vector  is  in  x-direction  Eq.  38  becomes 

1  dU  2 


,2  = 


_  fe  1  +  e  M  -  u 

2  dE  lcl  c2  J 

and  Eq.  40  yields: 


1  v  .2,  dU 
dE 


uv  =  j  (e^-  e22) 


(41) 


(42) 


Equation  41  and  42  are  the  well-known  equations  on  calculating  uv  and 
v2  when  the  velocity  vector  is  in  the  x-direction.  Put,  in  our  study 
when  the  wind  is  flowing  over  the  hill  the  velocity  vector  may  deviate 
from  the  x-direction.  Therefore,  Eqs.  38  and  40  were  1  sed  to  calculate 
the  uv  and  v2  when  $  ^  0  .  The  angle  4>  at  on  point  was 

estimated  from  the  streamline  pattern  shown  in  Fig.  5.  How  the  stream¬ 
lines  were  determined  will  be  discussed  later. 


For  crossed  wires,  when  <♦>  =  0  i.e.,  when  the  rlocity  vector  is 

parallel  to  the  x-axis,  the  angles  of  inclination  between  wire  1  and 
wire  2  and  the  x-axis  are  the  same  and  equal  to  45  (b  fih  wires  were  very 
carefully  mounted  perpendicular  to  each  other) .  In  or  jr  to  make  sure 
that  both  wire  1  and  wire  2  were  held  under  45°  to  x-a  Is  during  the 
experiment,  first,  the  wires  were  held  in  the  free  str  im,  when  the 
(ambient)  velocity  is  in  x-direction.  Then,  the  cross  1  wires  were 
rotated  180  about  the  hot  wire  probe  axis.  If  the  ou  nuts  of  the  wires 
were  different  after  this  rotation,  an  adjustment  in  t  :j  angle  of  the 
probe  with  flow  direction  was  made  until  the  anemometer  readings  of  both 
wires  were  invariant  to  rotation  about  the  longitudinal  axis. 


The  block  diagram  of  the  set  up  for  measuring  the  turbulence  is 
shown  in  Fig.  8.  A  single  wire  and  a  crossed  wire  were  mounted  side  by- 
side  at  the  same  height  on  a  24"  vertical  carriage.  The  elevation  of 
wires  could  be  read  off  as  a  voltage  across  a  potentiometer  geared  to 
the  positioning  shaft  and  was  either  read  out  from  a  digital  voltmeter 
(UVM)  or  plotted  on  an  x-y  plotter.  At  each  section,  data  were  taken 
at  55  test  points  shown  in  Fig.  9.  Also,  at  each  station  x  =  constant 
continuous  data  profile  plots  were  obtained  on  an  x-y  plotter.  The 
test  points  were  chosen  so  that  they  included: 

a.  points  on  the  trajectory,  i.e.,  points  on  the  trajectories 
at  the  distances  x  of  the  measuring  stations, 

b.  points  near  where  the  maximum  change  of  rms  value  of  u' 
occurred  in  each  section. 

At  each  of  the  test  points  the  following  data  were  taken: 

a.  The  rms  values,  i.e.,  the  fluctuations  in  voltage  of  a  single 
wire  and  of  two  crossed  wires.  All  three  rms  values  were  re¬ 
corded  by  x-y  plotters  versus  time  and  were  also  read  directly 
from  true  rms  meters  as  a  reference. 

b.  5-minute  turbulence  recordings  for  energy  spectrum  and  proba¬ 
bility  analysis.  A  Mincom  (Type  C100)  7  channel  FM  tape  recorder 
was  used  to  record  the  turbulence  for  both  single  and  crossed 
wires  (3-channel  simultaneous  recording).  The  output  of  the 
CSU-made  hot  wire  anemometer  has  a  dc  level  of  one  volt  and  an 
rms  value  of  the  order  of  0.05  volt.  The  dc  level  was  too  high 
and  the  rms  value  too  low  for  best  operation  of  the  tape  recorder. 
Therefore,  an  ac-amplifier  was  used  to  amplify  the  fluctuating 
voltage  and  to  eliminate  the  dc  level.  Furthermore,  an  atten¬ 
uator  was  connected  between  the  amplifier  and  the  tape  recorder 

to  adjust  the  recording  voltage  to  0.5  volt  rms. 

The  interconnections  of  all  instruments  are  shown  in  Fig.  8.  The  rms 
values  of  the  wires  obtained  from  the  rms -meter  (RMS  II)  before  the 
amplifier  and  attenuator  (A+A) .  The  recording  voltage  was  read  from 
RMS  II  of  Fig.  8. 

Besides  the  data  which  have  been  taken  at  each  test  point,  the  con¬ 
tinuous  rms  values  for  all  three  wires  were  also  recorded  on  an  x-y  plotter. 
Figure  10  is  a  typical  continuous  rms  profile  of  e  fur  a  single  wire 
at  x  =  12"  .  For  the  same  station  the  rms  profilis  of  e^  end  e^ 

for  wire  1  and  wire  2,  respectively,  are  shown  in  Figs.  II  and  12. 

As  long  as  the  rms  values  for  single  wires  and  corssed  wires  were 
known,  the  turbulence  components  ujv'  and  the  turbulent  shear  stress 
uv  could  be  calculated  from  Fq.  3,  £q.  11  and  Eq.  13,  respectively. 

The  measured  rms  values  and  the  calculations  of  u2  ,  v2  and 
uv  are  shown  in  Table  II.  The  profiles  of  these  quantities  are  plotted 
in  Fig.  13. 


21 


3.4  Determination  of  Streamline  Locations 


Mean  streamlines  can  be  drawn  so  as  to  be  always  tangential  to  the 
vectors  of  fluid  velocity  in  a  flow.  Since  a  separation  bubble  existed 
near  the  downstream  side  of  the  model,  it  was  found  desirable  to  first 
determine  a  reference  streamline  in  the  outer  part  of  the  flow  by  joining 
the  direction  of  mean  velocity  vectors  from  station  to  station,  and  to 
obtain  lower  streamlines  by  integration,  i.e.,  the  integral  \u  dy  =  constant 
below  the  reference  streanline  defined  other  streamliner.  ' 

The  first  streamline  was  found  by  using  a  hot  wire  in  the  following 
arrangement.  The  heat  transfer  from  hot  wire  to  air  depends  not  only  on 
the  magnitude  of  the  velocity,  but  also  on  the  flow  direction  with  respect 
to  the  wire.  The  heat  transfer  from  the  wire  is  maximum  when  the  flow  is 
perpendicular  to  the  wire,  minimum  when  the  wire  is  parallel  to  the  flow. 

By  rotating  the  hot  wire  and  plotting  the  output  of  the  hot  wire  anemometer 
versus  the  rotating  angle  on  an  x-y  plotter,  a  well-defined  minimum  was 
found  which  could  be  used  to  define  the  flow  direction. 

The  starting  point  of  the  reference  streamline  was  arbitrarily  set 
at  xq  =  0,  =  9.1"  ,  the  height  of  the  second  point  was  estimated  by 

y,  =  y  -l  sinAa  (43) 

1  1  J  o  o 

Aa  =  the  direction  of  local  velocity  vector  to  the 
free  stream  vector  at  first  point 

l  =  the  horizontal  distance  from  the  first  point  to 
the  second  point 

yQ  =  height  of  the  first  point 

In  general  the  height  of  the  point  yn  at  station  x^  can  be  cal¬ 
culated  from  the  n-lth  point  when  the  height  y n  ^  ,  the  angle  Ao^  ^ 

and  the  distance  i  ,  are  known.  We  have 

n- 1 


yn  =  Vi 


Z  ,  sinAa  , 
n-1  n-1 


(44) 


The  first  reference  streamline  was  estimated  in  this  manner,  with  the  re¬ 
sult  shown  on  Fig.  5. 

It  is  evident  that  this  is  not  too  satisfactory  a  method  for  deter¬ 
mining  the  streamline.  The  method  was  difficult  and  time  consuming,  and 
wrought  with  error  due  to  the  fact  that  the  error  in  estimating  the  height 
was  cumulative.  Therefore,  a  different  method  was  used,  where  integration 
starts  at  the  floor.  It  is  clear  that  outside  of  the  separation  region 
the  streamlines  can  be  determined  from  velocity  profiles  by  integrating 
up  from  the  lower  boundary. 


According  to  continuity  the  flow  rate  between  two  streamlines  at  any 
section  must  be  the  same.  It  is  known  that  the  lower  boundary  is  a  streamline. 
The  lower  boundary  consists  of  three  parts 

(a)  before  separating  the  surface  of  the  model  hill,  upstream 
of  separation 

(b)  between  separation  and  reattachment,  the  upper  boundary 
of  the  separation,  bubble 

(c)  after  reattachment,  the  floor. 

In  our  flow  field  the  part  (a)  and  (c)  are  fixed  and  are  well  known.  Thus, 
a  reference  streamline  can  be  found  for  regions  (a.)  and  (c)  at  some  height 
up,  and  the  streamline  above  the  separation  region  (which  is  rather  short) 
can  be  found  by  fairing  between  the  l»o  curves  following  the  trend  of 
the  reference  streamline  measured  by  the  previous  method.  The  remainder 
of  the  procedure  of  stream.’  ne  construction  is  as  was  outlined  above. 

In  this  manner  two  streamlines  are  obtained,  namely,  (1)  the  measured 
streamline  starting  5"  above  the  crest,  which  was  constructed  from  the 
reference  str'eamline  by  applying  the  continuity  equation,  and  (2)  the 
lower  boundary.  Between  these  two  streamlines  the  flow  rates  are  q 
at.  any  section.  Some  streamlines  were  interpolated  between  these  two. 

The  streamlines  and  the  separation  region  are  shown  in  Fig.  5. 

3.5  Measurement  of  Turbulence  Spectra 

Spectra  of  the  u' -component  of  the  turbulent  signal  was  obtained 
bv  means  of  a  Bruel  and  kjaer  spectrum  analyzer  (Type  B  5  K  2109),  with 
occasional  cross  checks  against  results  from  a  Technical  Products  Wave 
Analyzer  (Type  TP  62).  The  former  has  a  proportional  band  width, 
passive  filter  system,  while  the  latter  works  with  a~*ive  constant 
banuwidth  filters.  Both  set-ups  for  this  evaluation  ,re  shown  in  Fig.  14. 

3.6  Measurement  of  Probability  Uenbi ties 

Ihe  probability  density  distribution  of  a  single  turbulent  component 
was  measured  with  a  Technical  Products  probability  analyzer  (Type  TP  647), 

Fig.  15.  Joint  probability  densities  were  measured  with  two  of  the  above 
analyzers  coupled  tosether  so  that  one  provided  the  gate  for  the  joint 
probability  density  obtained  from  the  other.  Fig.  16. 

Normalization  and  calibration  of  probability  analyzers  were  based 
on  a  known  input  sine  wave,  whose  rms  value  is  close  to  that  of  the  hot¬ 
wire  turbulent  signals.  The  probability  density  or  a  sine  wave  e=AsinO 
whose  phase  0  is  a  random  variable  uniformly  distributed  on  the  interval 

i!  11  .  .  . 

-T  to  T  ,  is  given  by: 


f  t  e ) 


for  e  -  A 


» 


0 


otherwise 


(45) 


Since  the  normalising  process  has  mads  all  the  amplitudes  of  different 
sine  waves  to  be  A  =  -Jl  Jo*-  =  /f  in  such  a  way  that  =  1 


the  lowest  point  of  the  sine-wave  probability  density  is  found  to  be 
f(o)  =  — ^-  =  0.225  ,  which  was  used  for  calibrating  the  x-y  plotter. 

7T  A 

An  example  is  given  in  Fig.  17.  As  for  the  calculation  after  the  analog 
analysis,  the  main  problem  was  to  convert  the  measured  a  values  into 
the  real  turbulent  fluctuation  units  (in  feet  per  second).  A  graphical 
integration  based  on  the  second-moment  of  the  probability  density  was 
suggested  as  a  proper  approach,  i.e,. 


u.  i 

l 


(46) 


where  u^':  is  the  square  of  the  rms  value  at  point  i  and  is  the 

value  obtained  from  the  probability  analyzer. 

The  experimental  data  for  probability  densities  were  plotted,  both 
for  probability  densities  of  single  quantities  and  for  joint  probability 
densities.  Instead  of  establishing  a  3-dimensional  distribution  of  joint 
probability  density,  iso-probabilitv  density  contour  maps  were  plotted. 
Conditional  density  functions  were  evaluated  according  to  the  definition, 


f(v'/u') 


f(u’,v') 


f 

i 


f (u* ,v')dv( 


(47) 


and  thus,  the  conditional  probability  density  is  given  as 


f(v'/u'  =  u^jdv'  - 


r 

Prob  u*  <  u'  <  u'  +  du'  v'<v'<v*  +  dv' 

L  o  _ ~  o  *  o  ~  o _ 

Prob i  u  <  u’  <  u'  +  du' 

(_  0  o 


(48) 


This  equals  the  ratio  of  the  mass  in  the  differential  element  of  Fig.  18 
to  the  mass  in  the  strip  (u^  ,  u^  +  du')  .  Thus,  for  a  given  , 

the  density  f(v'/u'  =  u^)  is  the  u^  -  profile  of  joint  density  f(u'v') 

normalized  to  make  its  area  equal  to  1  ,  (Fig,  18). 

3 . 7  Measurement  of  Space  Correlation  Coefficients  Along  the  Trajectories 

Space  correlations  along  trajectories  were  taken  by  passing  the 
outputs  of  two  single  hot-wire  anemometers  through  a  well-calibrated  sum- 
and-differenc.e  circuit  instead  of  an  analog  multiplier  (Fig.  19).  The 
calculations  were  only  based  on  the  rms  values  of  inputs  and  outputs  of  the 
sum-and-di fference  circuit. 
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where  As  ,  A^  ,  are  the  calibration  constants  for  hot-wire  anemometers, 
1  and  2,  respectively,  Cg  ,  ,  the  calibration  constants  for  sum  lS) 

and  difference  (D)  circuits.  Then 
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yields  a  relation  for  the  space  correlations  when  £  is  the  distance  of 
point  2  from  point  1,  and  x  is  the  location  of  point  1  in  the  reference 
coordinate  system. 
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Chapter  IV 


THE  EXPER CENTAL  RESULTS 


The  experimental  work  on  this  project  was  conducted  in  three  I 

phases.  These  were:  I 

i 

1.  Measurements  of  the  mean  wind  vertical  velocity  profiles  and 
turbulenf  intensities  at  selected  points  on  the  lee  side  of  a  sinusoidal 
hill  using  the  Army  Wind  Tunnel.  This  work  has  been  reported  by  Plate 
and  Lin  (1965),  .Another  more  detailed  flow  pattern  will  be  given  in 
this  chapter. 

2.  Determination  of  theoretics1  missile  trajectories,  corresponding 
to  mean  wind  conditions  if  the  missiles  were  fired  from  the  lee  side  of  a 
scaled-up  version  of  the  two-dimensional  ridge. 

3.  Determination  in  the  wind  tunnel  of  the  characteristics  of 
the  wind  field  at  selected  points  along  the  scaled-down  missile  tra¬ 
jectories,  in  accordance  with  the  theoretical  development  given  in 
Chapter  II. 

A  fourth  phase,  not  reported  here,  will  be  to  calculate  the  response 
of  the  missile  to  the  experimental  wind  fieMs  determined  in  phase  3. 

4 . 1  Determination  of  Mean  Missile  Trajectories 

This  work  was  conducted  at  the  USA  Ballistic  Research  Laboratories, 

Aberdeen,  Maryland,  using  the  laboratory  compiling  facilities  and  a  six 
degree  of  freedom  multi-stage  rocket  trajo  t.ory  program. 

The  missile  used  for  the  mean  trajectory  calculation  was  a  hypo- 
t'  etical  gun-launched  two-stage  anw-tank  missile.  The  gun  launched 
the  missile  at  1200  f/s.  After  a  short  delay,  a  booster  section  ignited, 
the  thrust  from  which  accelerated  the  missile  to  a  velocity  of  2100  f/s. 

At  that  point,  a  sustainer  motor  ignited,  the  thrust  from  which  kept  the 
missile  at  a  constant  velocity  until  it  re ache  1  a  position  about  1  km 
from  the  launcher.  In  computer  simulations,  this  missile  was  shown  to 
have  a  steady  cross  wind  sensitivity  of  0.36  mils  angular  deflection  per 
ft/sec  of  cross  wind. 

For  the  simulation  study,  the  twodi-.nsional  ridge  used  in  the 
tunnel  was  scaled  up  by  a  factor  of  1200  to  a  ridge  400  ft  high  by  2000  ft 
long.  It  was  then  assumed  that  missile  launchers  were  emplaced  at  the 
base  of  the  ridge;  halfway  up  the  ridge,  and  at  the  top  of  the  ridge. 

All  the  launchers  were  pointed  at  targets  oi.  the  lee  side  of  the  ridge, 
the  targets  being  1  km  from  the  launcher  sites. 
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The  trajectories  of  missiles  were  simulated  first  for  the  no-wind 
case  and  then  for  the  case  of  the  steady  wind  flowing  over  the  ridge  by 
interpolating  in  the  data  from  Plate  and  Lin  (1965) . 

The  missile  trajectory  data  from  these  simulations  were  then  sent 
to  Colorado  State  University  to  be  used  in  further  experimental  work. 
The  characteristics  of  the  wind  fields  along  the  trajectories,  shown  in 
Pig-  3,  generated  in  the  above  study,  is  discussed  in  this  report. 

4.2  Me^.  Velocities  and  Streamline  Pattern 


Mean  velocity  distributions  are  shown  in  Fig,  5.  The  solid  lines 
indicate  velocities  measured  with  a  pitot  static  tube,  while  the  dashed 
lines  refer  to  hot-wire  measurements.  On  the  whole,  the  agreement  be¬ 
tween  the  two  sets  is  good,  even  without  any  corrections  for  turbulence. 

The  small  deviations  might  just  as  well  be  due  to  drift  in  the  hot-wire 
characteristics,  which  could  never  be  fully  eliminated. 

Characteristic  of  the  flow  field  is  the  strongly  accelerated  riuw 
above  the  crest  of  the  model,  which  gives  rise  to  the  velocity  maximum, 
and  the  very  sharp  velocity  gradients  in  the  neighborhood  of  the  separation 
streamline.  These  velocity  gradients  interact  with  the  turbulent  shear 
stress  to  cause  a  large  increase  in  the  amount  of  turbulent  energy  of 
the  f  1  -w . 

Vertical  mean  velocities  can  be  determined  from  the  mean  streamline 
pattern  shown  in  Fig.  5.  The  streamline  pattern  also  shows  the  separation 
region.  Under  the  separation  streamline,  the  velocity  gradually  decreases, 
reaches  zero  and  reverses  direction.  This  can  be  inferred  from  the  fact 
that  the  discharge  across  any  vertical  section  underneath  the  separation 
streamline  must  be  zero.  The  experimental  data,  however,  fail  to  show 
this  behavior  due  to  the  fact  that  the  pitot  tube  cannot  measure  any 
backflows,  while  the  hot-wire  cannot  distinguish  the  flow  directions. 

4.3  Turbulent  Intensities  and  Shear  Stresses 

The  turbulent  quantities  u'2  and  v'2  ,  and  the  turbulent  shearing 

stress  u*v'  were  plotted  against  y  in  Fig.  13.  The  turbulent  quantity 
w'2  was  also  determined  at  a  number  of  points  and  generally  behaved  roughly 
like  the  u-component.  The  shear  stresses  u  w  were  calculated  and  found 
to  be  negligibly  small  even  in  the  bubble  region.  The  profiles  shown  in 
Fig.  13  have  a  strongly  peaked  shape  in  the  neighborhood  of  the  separation 
streamline,  especially  for  short  distances  from  the  separation  point  on 
the  hill  slope.  In  general,  the  intensity  profiles  of  u  component  based 
on  values  which  had  not  been  corrected  for  the  flow  direction  coincide 
with  those  obtained  by  Plate  and  Lin  (1965)  and  thus  the  reproduct ibility 
of  the  turbulent  flow  field  in  the  present  tunnel  is  quite  satisfactory. 

4 . 4  Turbulence  Spectra  and  Dissipation  Rates 

Turbulence  spectra  were  evaluated  for  all  points  indicated  in  Fig.  9. 
The  signal  is  plotted  in  the  form  eB2/B  vs.  f  (Fig.  20  through  Fig.  28). 

Here  eQ2/B  is  the  energy  density,  per  Hz  ,  of  the  electrical  signal 
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from  the  hot  wire  as  passed  through  the  filter  of  bandwidth  B  and  of 
center  frequency  f  of  the  spectrum  analyzer.  It  differs  from  the 
energy  level  of  the  turbulent  motion  by  a  calibration  factor  given  by 
the  square  of  the  slope  of  the  calibration  curve  U  vs.  E  of  the  hot¬ 
wire  anemometer. 

In  the  low  frequency  range,  we  notice  a  strikingly  different 
spectrum  shape  close  to  the  hill  crest  (x  =  0)  as  compared  to  the 
results  at  16  inches  downstream.  At  short  vertical  distances  from  the 
wall,  the  data  close  to  the  crest  (Fig,  22)  indicate  a  much  slower  drop¬ 
off  with  frequency  than  the  set  of  data  shown  in  Fig.  25.  In  fact,  there 
seems  to  exist  a  well-developed  region,  between  40  and  200  cps,  in  which 
the  energy  level  decreases  almost  linearly.  This  behavior  is  character¬ 
istic  of  strong  interactions  between  mean  flow  and  turbulence,  i.e., 
of  a  flow  when  a  large  amount  of  turbulence  generation  due  to  large 
velocity  gradient  takes  place.  This  behavior  is  not  typical  for  other 
boundary  layer  flows  of  the  U.S.  Army  Wind  Tunnel. 

Due  to  strong  noise  levels  of  the  magnetic  tape  recorders,  the 
part  of  the  spectra  corresponding  to  frequencies  above  2000  Hz  is  not 
usable.  For  large  frequencies,  but  below  2000  Hz,  the  shape  of  the 
spectrum  is  the  same  for  all  data.  In  fact,  if  the  spectrum  is  plotted 
in  the  similarity  form  of  the  universal  equilibrium  law  of  Kolmogoroff , 
we  find  that  the  shape  is  identical  for  all  data,  and  they  collapse  on 
a  single  curve.  This  is  illustrated  in  Figs.  29,  30,  31,  32  and  33,  in 
which  the  data  of  Figs  22,  23,  25,  26  and  28  have  been  replotted  in 
dimensionless  form: 

*(£-)=  N.F(f)  (50) 

s 

<p  is  the  non-dimensional  spectral  density,  F(f)  is  the  measured 
spectral  density  at  frequency  f  ,  and  N  is  a  conversion  factor: 
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Furthermore,  K  is  the  reference  wave  number  based  on  the  dissipation 
c  : 


..  ,  -3  1/4 

KS  =  (rv  } 


(52) 


As  an  estimate  for  the  dissipation  e  we  have  used  the  isotropic 
relationship: 

I  ")U '  1  ?  -2  3u '  1  2 

f  -  15  “  iirj  ■ 15 » "  irj  (53> 

and  also,  the  equivalent  form 
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where  K  is  the  wave  number 


K 


2  TT  f 
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(55) 


and  D"  is  the  mean  velocity,  as  before.  Dissipation  rates  calculated 
from  Eqs.  (53)  and  (54)  are  shown  in  Fig.  34. 

In  Figs.  29 t  30  f  31(  32  and  33  we  have  also  indicated  the  -5/3 
law  of  the  inertial  subrange  and  the  universal  shape  of  the  high  fre¬ 
quency  and  of  the  turbulence  spectrum,  in  the  form  given  by  Sandbom 
and  Marshall  (1963).  It  is -surprising  to  see  that  the  high  frequency 
end  of  the  -spectrum  in  the  highlv  disturbed  boundary  layer  of  our  case 
is  presented  exactly  by  the  high  frequency  shape  of  the  undisturbed 
turbulence  in  a  boundary  layer  along  a  flat  plate.  Since  Sandbom  and 
Marshall  have  demonstrated  the  perfect  agreement  of  their  spectra  with 
experimental  results  obtained  in  wind  over  ocean  waves  by  Pond  et  al., 
(1963),  it  appears  that  this  range  of  the  spectrum  is  a  universal  feature 
of  all  turbulent  flows. 


But  the  same  conclusion  cannot  be  drawn  for  the  turbulence  spectrum 
in  the  inertial  subrange.  Pond  et  al.,  report  that  here  a  spectrum  law 
of  the  form  is  valid  ' 


N  F(f)  =  k  (£-)  (56) 

s 

where  k  is  a  universal  constant,  about  0.46.  For  the  data  of  this  re¬ 
port,  it  is  found  very  near  the  crest  of  the  model  hill,  this  "constant" 
is  well  enough  verified,  but  at  larger  distances  downstream,  in  the 
region  which  derives  its  turbulence  from  the  initially  strong  gradients 
in  the  mean  velocity  across  the  separation  streamline,  the  "constant" 
seems  to  be  substantially  higher.  At  16"  (Fig.  31)  downstream  from  the 
hill  ''rest,  the  best  fitting  -5/3  law  has  a  constant  k  of  about  0.85. 
It  should  be  noted  that  in  the  velocity  region  where  this  is  found  the 
turbulence  level  decreases  rapidly  with  distance,  indicating  that  the 
amount  of  energy  generated  locally  is  lower  than  that  dissipated,  i.e., 
the  ratio  of  dissipation  to  generation 


D 


(57) 


in  this  region  is  greater  than  one.  Inis  result  thus  is  in  qualitative 
agreement  with  a  result  of  Margolis  and  Lumley  (1964).  It  has  as  yet, 
however,  not  been  shown  that  3  universal  relation  exists  between  k  and 
D  .  Experiments  are  at  present  underway  at  Colorado  State  University 
to  investigate  this  point.  That  D  might  also  be  an  important  quantity 
in  modeling  of  atmospheric  turbulence  has  been  pointed  out  by  Plate  and 
Lin  (1966). 
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The  low  frequency  end  of  the  spectrum  is  governed  by  the  process 
of  energy  extraction  from  the  mean  flow  and  depends  on  the  local  velocity 
field.  Similarity  forms  can,  therefore,  not  be  expected  for  the  whole 
spectrum.  But  the  eddies  associated  with  the  low  frequency  end  of  the 
spectrum  input  cause  the  most  important  dynamic  effect  on  a  missile  during 
its  flight.  Work  is  therefore  in  progress  at  CSU  on  relations  between  the 
low  frequency  end  of  the  spectrum  and  the  local  mean  velocity  field. 


4 . 5  Probability  Density  Pistribut ions 

Two  different  sets  of  probability  densities  are  given  for  the 
points  along  the  trajectories.  These  are  probability  densities  of  the 
v',  of  the  ’  components  (Figs.  35-42),  and  joint  probability  densities 
of  u'  and  v'  , (Figs.  43-50). 

The  density  of  each  individual  turbulent  component  seems  quite 
well  represented  bv  a  Gaussian  (normal)  curve  except  the  mean  value  is 
not  exactly  zero.  This  is  illustrated  by  plotting  cumulative  probability 
densities  on  probability  papers  (Figs.  51-54).  The  angle  9  of  the  joint 
probability  density  contour  as  calculated  by  Eq.  (26),  which  is  based  on 
the  assumption  of  a  zero  mean  density,  is  surprisingly  close  to  the 
measured  angle  6  obtained  directly  from  the  two-dimensional  contour  of 
joint  probability  density  (Table  III).  Most  of  probability  densities 
also  show  the  evidence  of  skewness.  Numerical  evaluation  of  the  skewness 
factor  is  not  necessary  for  this  study,  but  the  reason  for  this  skewness 
could  always  be  interpreted  as  the  result  of  the  preferred  direction  of 
fluctuation  of  the  turbulence  component  due  to  the  shear  stresses.  This 
phenomenon  could  be  strongly  show'’  for  those  points  along  the  separation 
streamline  and  somewhere  within  the  bubble  region  because  of  the  complex 
interactions  of  the  three-dimensional  turbulence. 

Flatness  of  the  probability  density,  which  was  suspected  as  the 
result  of  turbulent  intermittency ,  was  not  obviously  seen  among  our 
measurements  which  were  usually  performed  in  the  fully  turbulent  boundary 
1 aye  r . 


In  order  to  check  the  statistical  independence  of  w'  from  u’  , 
joint  probability  density  distributions  f(u',  w')  are  measured.  As 
mentioned  in  a  previous  section,  the  homogenuity  of  turbulence  in 
plane  parallel  to  the  ground  squires  that  u'w’  be  zero,  which  was 
confirmed  by  our  experiments.  But  it  is  also  necessary  that  for  each 
value  of  u'  ,  f(w''u')  be  independent  of  u'  and  equal  to  f(w'). 

To  check  this  criterion  (which  gives  a  sufficient  condition  for  sta¬ 
tistical  independence)  we  have  plotted  into  one  graph  f(w'|u')  for 
different  values  of  u'  as  well  as  f(w')  .  Some  results  are  shown  in 
Figs.  55-59.  They  indicate  that  u'  and  w'  are  not  exactly  jointly 
Gaussian,  even  though  u1  w'  was  found  to  be  zero  from  hot  wire  measure¬ 
ments,  which  is  not  a  sufficient  condition  to  conclude  that  w'  is  inde¬ 
pendent  of  u'  for  the  slightly  non-Gaussian  curves  observed  in  this 
study . 


We  next  checked  joint  probability  density  distributions  f(u'  u'  ) 

'  ii-l 

different  results  were  obtained  ^upending  on  the  direction  of  the  tra¬ 
jectory.  For  a  trajectory  of  60  azimuth  these  typical  measurements 
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are  shown  in  Figs.  60-62.  The  results  indicate  that  no  correlation  exists 
for  turbulence  components  u',  and  u'^  .  This  is  evident  from  the  co¬ 

incidence  of  the  axes  of  the  elliptical  contours  with  one  of  the  axes.  It 
implies  that  the  distance  between  the  two  points  has  been  chosen  too  wide, 
so  that  the  integral  length  scale  is  somewhat  smaller  than  the  distance 
between  adjacent  points. 

For  the  case  of  60°  azimuth,  the  assumption  of  an  elongated  eddy 
structure  with  u^ '  related  only  to  u., '  is  not  justified.  The  only 

way  of  obtaining  an  impact  dispersion  probability  distribution  is  found 
by  assuming  statistical  independence  between  neighboring  points  along 
the  trajectory,  with  the  resulting  inconsistency  in  the  averaging  pro¬ 
cedure  for  the  average  velocity  encountered  by  the  missile  in  traveling 

from  x.  ,  to  x. 

l-i  i 

Whan  joint  probability  densities  for  ul  ^  and  u.  1  lying  on  a 

trajectory  in  the  mean  wind  direction  (0°  azimuth)  are  considered,  strong 

correlations  are  found  as  is  demonstrated  in  Figs.  63-65.  The  results 

are  even  better  illustrated  by  considering  directly  the  space  correlation 

functions  along  different  trajectories,  which  are  plotted  in  Fig.  66. 

Two  trajectories  were  added,  trajectory  2  which  starts  halfway  down  the 

hill  slope,  and  trajectory  3  which  starts  at  the  crest  of  *he  ridge. 

For  correlation  functions  which  start  at  the  hill  directly,  correlation 

distances  are  rather  short:  ai ready  at  1"  downstream  from  the  hill  the 

correlation  between  u  and  u,„  has  decreased  to  a  correlation  co- 

o  1" 

efficient  r  =  0.2  ,  where 


For  correlations  starting  at  larger  values  than  x  =  4"  the  correlation 
coefficients  for  both  trajectories  seem  to  follow  a  fairly  constant 
correlation  curve,  with  approximately  equal  shape  and  furtnermore,  with 
a  shape  which  is  approximately  exponential.  The  exponential  curve  which 
appears  to  give  the  best  fit  has  been  sketched  into  the  experimental 
results.  It  is  represented  by  the  equation: 

r  a  e"aX  (59 ) 

when  a  is  a  constant  equal  to  0.2  in  1 .  This  result  implies  that  in 

the  outer  region  (further  downstream  than  4"  from  the  ridge)  the  joint 

probability  distribution  f  ( u  ? ,  u.'  )  can  be  specified  by  the  magnitudes 

7  7  11-1 

u!‘  and  u!“  and  by  the  curve  expressed  by  Eq.  (.39:  If  u'  and  u' 

can  be  approximately  represented  by  Gaussiani ty .  The  latter  yields  the 
numerator,  and  the  former  the  denominator  of  f.q .  (26)  with  u'  ^  taking 

the  place  of  v  , while  u'  takes  the  place  ot  u1 .  \  difficulty  in 

interpretation  arises,  however,  if  u.  happens  to  be  equal,  or 
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nearly  equal  to  u.,  ^  .  In  that  case,  we  may  let  u+  =  Au^  and 

u+  =  j  ,  where  A  is  any  positive  constant  factor  in  order  to 
either  enlarge  or  contract  the  values  of  in  such  a  way  that  u+ 

and  u+  j  are  numerically  quite  different  from  Eq.  (21)  and  assuming 
again  m  =  m  =0  for  simplicity, 


u. 

1 

Vi 

f(u+ 

1 

*  ut-l 

ftv  VP 
J  <v  vP 


.  u+ 

=  I  H7T  ’  ut-l5 


2ttAo  o  i 
u.  u.  , 

1  1-1 


'1-p2 


t 

1 


2(l-p2) 


u+ 

1 


A  a 


L* 


-  2P 


i-1/  . 


>(60) 


Now,  following  the  same  procedure  in  deriving  our  Eq.  (26),  we  may  obtain 


Tan  20 


2A  u. 

l 


u. 

l- 


1 


u2 


i-1 


(61) 


But  this  angle  0  can  only  be  used  to  specify  f(u+,  u+  p  ,  and  of 
course,  f(u^,  u^  p  can  also  be  specified  by  a  little  transformation 
technique . 

For  the  region  between  0  and  4"  along  the  trajectory,  it  is 
necessary  to  know  the  coefficient  a  for  each  length  section,  which  can 
approximately  be  found  by  linear  interpolation.  It  is  thereby  assumed 
that  the  correlation  function  is  exponential  between  point  x^  and 

x^  1  .  Thus,  even  though  at  larger  distances  the  same  exponential  does 

not  fit  all  the  data,  it  is  possible  to  make  the  assumption  of  Markoff 
dependency  over  a  short  reach,  with  at  least  engineering  accuracy. 
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Cone lusions 


In  the  report  we  have  given  experimental  evidence  that  the 
turbulence  structure  in  the  complicated  flow  field  downstream  of  a 
model  hill  can  be  represented  to  a  good  approximation  by  Gaussian  and 
jointly  Gaussian  probability  distributions  connected  by  space  correla¬ 
tion  functions  which  appear  exponential.  The  joint  proba  ..  ity  den¬ 
sities  are  seen  to  be  fully  specified  by  the  variances  of  the  two 
variables  ani  their  first  order  cross  correlation. 


The  vesults  of  this 
construct  a  model  of  the  • 
by  simply  determining  the 
towards  this  objective  is 


-udv  indicate  that  it  is  possible  to 
.ssi le  dispersion  probability  distribution 
variance  of  this  distribution.  Research 
currently  under  wav  at  CSU. 


!< 
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TABLE  I 


Mean  Velocity  Calculation 


p* 

y 

in 

!J 

HW 

UP.T 

X 

in 

pt 

y 

in 

UHW 

UP.T 

6 

3.89 

22.4 

29.4 

8" 

18 

1.35 

3.6 

7 

3.98 

29.0 

27.8 

19 

1.52 

4.4 

~ 

8 

4.09 
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33.3 

20 

1  .9 

7.9 

10.1 

9 

4.22 

32.0 

33.3 

21 

2.0 

10.6 

1 1 . 9 

10 

4.68 

31.8 

32.9 

22 

2.16 

13.6 

15.5 

11 

5.1 

32.0 

25 

2.51 

19.7 

20.2 

6.0 

31.4 

3  x ,  3 

24 

3.12 

22.6 

23. 3 

8.0 

31.2 

30.8 

25 

4.86 

24.7 

25,3 

10  0 

:>i . 0 

30.7 

6. 

26.3 

27.4 

15.0 

31.0 
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8. 

27,2 

28.3 

20.0 

30  9 

30.6 

10. 

28.5 

29.3 

24.0 

30.7 

30.2 

15. 

29.5 

30.2 

20. 

30.3 

30.6 

12 

3.35 

16.2 

19.2 

24. 

30.3 

30.4 

13 

3.49 

26.5 

?C  1 

14 

3 . 60 

27.3 

D 

10" 

0.05 

1.15 

1  .707 

IS 

3.84 

29.0 

29.8 

.7 

5.15 
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16 

4.63 

29.0 

29.7 

1.15 

7.96 

8.18 

17 

5.25 
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29.7 

1.65 

16.8 

12.77 

6.0 

29,5 

30 . 1 

2.2 

22.5 

17. 

8.0 
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3 . 36 
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30.5 
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5.43 

25.8 
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15.0 

30 . 5 
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8. 
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26.72 
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31. 

30.5 

10. 

28.5 

28.2 

24.0 

30.7 

30.5 

15. 

29.5 

30.06 

20. 

30.6 

31. 

2  46 

10.2 

2. S3 

16.6 

8.18 

12" 

26 

0.12 

3.8 

- 

2.81 

24.1 

1 7 .  u6 

27 

0.73 

5.9 

5.12 

3.29 

26.2 

23.22 

28 

1.12 

8.0 

9.2 

4.0 

27.4 

24.58 

29 

1.30 

1 1  .1 

13.3 

6. 

28.6 

26.77 

30 

1.75 

14.4 

15. 

8. 

29.7 

28,26 

31 

2.66 

20.0 

20.9 

10. 

30.4 

29.2 

32 

5.00 

24.0 

26.2 

15. 

51 .0 

30.5 

7.14 

25.4 

25.6 

20. 

31  .4 

31.16 

9. 

27. 

27.4 

12. 

28.2 

28.5 

15. 

28.9 

29 . 5 

23.3 

29.9 

30.3 
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X 

in 

Pt 

y 

in 

UHW 

llP.T 

14” 

0.05 

2.7 

.18 

4.14 

- 

.4 

6.0 

6.15 

.62 

7.96 

8.86 

1.05 

11 . 

12.31 

1.58 

16,6 

15.55 

2.0 

19. 

17.16 

2.67 

21.35 

19.62 

4. 

22.64 

22.44 

8. 

26.8 

26.4 

10. 

27.5 

27.64 

iS. 

28.9 

29.58 

20. 

30.4 

30.64 

16” 

53 

0.12 

5.8 

9.66 

74 

.28 

8.4 

11.32 

35 

.69 

10.2 

13.87 

36 

.96 

.10.8 

15.17 

37 

1  , 56 

13.8 

18.07 

38 

2.2 

18.2 

19.62 

39 

3.52 

20. 

22.6 

40 

6.5 

23.8 

25.9 

8. 

25.8 

26.98 

10. 

27,8 

28.5 

15. 

28.8 

30. 

20. 

29.9 

30.7 

18” 

.05 

7.4 

9 , 66 

.13 

9.55 

11.45 

.37 

10.95 

12.66 

.95 

14.05 

13.44 

2. 

19.5 

17.4 

2.9 

20.85 

20.42 

4. 

22.3 

22.6 

8. 

26.3 

26.42 

10. 

27.7 

27.76 

15. 

29.5 

29.6 

20. 

30 . 2 

30 . 64 

an  tin  lied) 


A 

.in 

pt 

3 

in 

UHW 

.T 

24” 

41 

0.12 

12. 

12.81 

44 

1.06 

15.8 

17.31 

45 

1.20 

16.2 

16.85 

46 

2.73 

18,8 

20.25 

47 

6.16 

23.9 

23.85 

8. 

25.7 

25.2 

10. 

26.4 

26.24 

12.5 

27.55 

27.65 

15. 

28.20 

28,58 

20. 

29.3 

29.4 

32” 

48 

0.14 

14. 

14.02 

49 

.50 

16.5 

16.85 

50 

3,58 

20.9 

21,82 

51 

6.74 

24. 

25 , 55 

10. 

26.55 

26.8 

12.5 

27.1 

27.3 

15. 

28.1 

28.56 

20. 

29.4 

29.48 

40" 

52 
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14.5 

15.4 

53 

.2 
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16.57 

54 

4.16 
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55 

6,32 
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25.1 
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Turbulent  Calculation 
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TABLL  II  (Continued) 


x  pt _ v _ u  ' 2 


6 

3.89 

27.7 

7 

3.98 

26.5 

8 

4.09 

3.37 

9 

4.22 

2.46 

10 
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1.88 

11 
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13 
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19 
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25 
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u' v' 

.227 

18.51 
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TABLE  II  (Continued) 
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6  =  1/2  tan 
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